Limits on Neutrino Mass from Cosmic Structure Formation 
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We consider the effect of three species of neutrinos with nearly degenerate mass on the cosmic 
structure formation in a low matter-density universe within a hierarchical clustering scenario with 
a fiat initial perturbation spectrum. The matching condition for fluctuation powers at the COBE 
scale and at the cluster scale leads to a strong upper limit on neutrino mass. For a flat universe 
with matter density parameter = 0.3, we obtain < 0.6 eV for the Hubble constant Hq < 80 
km s~^ Mpc~^. Allowing for the more generous parameter space limited by f2 < 0.4, Ho < 80 km 
s^^ Mpc^^ and age to > 11.5 Gyr, the limit is 0.9 eV. 
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Recent experiments for atmospheric and solar neutrino 
fluxes suggest that the neutrinos are massive. In partic- 
ular, the atmospheric neutrino experiment indicates an 
ahuost maximal mixing between the two neutrinos, which 
is most naturally understood if the relevant species are 
nearly degenerate in mass. Nearly maximal mixing is 
also a viable possibility to explain the long-standing so- 
lar neutrino problem with oscillation either in vacuum or 
in matter, although there remains the solution that it is 
explained by small-angle mixing via oscillation in matter 
[1] . For these reasons the idea has gained popularity that 
the three neutrinos are massive and almost degenerate in 
mass (e.g., [1,2]). The degenerate neutrinos mean that 
neutrino mass is larger than several tenths of eV, and 
this means that they provide the universe with a matter 
density comparable to or more than that in stars, and 
play some role in cosmological structure formation. 

There are a few authors who discussed the possibility 
that neutrinos have played an active role in the forma- 
tion of large-scale structure of the universe, especially in 
giving a power at a large scale which otherwise cannot 
be accounted for in the standard cold dark matter sce- 
nario at the critical matter density [3] . At the time of the 
emergence of this idea theorists took more seriously the 
Einstein-de Sitter (EdS) universe of the critical matter 
density, so that typical compositions of the matter were 
assumed to be VLcdm = 0.7 — 0.8 and Vl^, — 0.3 — 0.2 in 
units of the closure density, 10.54ft,^keV (cm)^'^, where h 
is the Hubble constant in units of 100 km s~^ Mpc~^. 
This neutrino mass density corresponds to neutrino mass 
of (30 — 20) ft.^ eV. There have been many explorations of 
this scenario since the proposal [4] , and the current con- 
clusion is that the neutrino density in excess of fi^ > 0.3 
is disfavoured in the EdS universe from the viewpoint of 
early cosmic structure formation. 

Over the last few years the evidence has been accu- 
mulated indicating a low density universe. There are 
also observations pointing to the dominance of the vac- 
uum energy (cosmological constant. A) that makes uni- 
verse's curvature flat, which is also preferred from the- 



oretical point of view for a low matter density universe. 
The list in favour of a low matter-density universe in- 
cludes; (1) Hubble constant - cosmic age mismatch for 
the Vl — \ universe; (2) No positive indications for the 
presence of copious matter beyond the cluster scale: the 
mass to light ratio inferred from clusters and galaxies, 
M/L = (100 - 400)/i, corresponds to Vt = 0.1 - 0.3 
[5]; (3) consistency of the cluster baryon fraction with 
the field value [6]; (4) The slow evolution of the cluster 
abundance from redshift z—Q to 0.8 together with the 
abundance normalization at z=0 [7]; (5) The Hubble di- 
agram of Type la supernovae [8]. The data indicate a 
finite cosmological constant. In view of systematic er- 
rors in various steps of the analysis, however, a zero A 
is probably not excluded, whereas an O = 1 universe is 
too far away from the observations; (6) The perturbation 
spectral-shape parameter F VLh — 0.2 — 0.3 from large- 
scale structure [9]; (7) Matching of the power spectrum 
between COBE and galaxy clustering [10]; (8) Evolution 
of small scale non-linear galaxy clustering [11]; (9) Lo- 
cal velocity field versus density enhancement [12]. The 
results of (1) — (9) converge between fio — 0.2 and 0.4. 

The two positive indications for the presence of a cos- 
mological constant are the Type la supernova Hubble 
diagram mentioned above and the acoustic peak distri- 
bution in the cosmic microwave background radiation 
(CBR) anisotropics [13], which says that the universe 
is close to fiat whether dominated by matter or vacuum 
energy. 

We remark that the results from large-scale velocity 
flow analyses are controversial; the resulting SI varies 
from analyses to analysis compounded by the uncertainty 
in the biasing parameter regarding the extent to which 
galaxies trace the mass distribution [14]. Recent analyses 
of galaxy peculiar velocities combined with other obser- 
vations claim that they are consistent with a low density 
universe with a finite A [15]. 

There seems no extensive analysis available for the ef- 
fect of neutrinos in a low matter-density universe, al- 
though a brief reference has been made in [4]. In this 
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paper wc consider the effect of massive, three degener- 
ate neutrinos on the cosmic structure formation in the 
low matter-density universe. We assume the hierarchi- 
cal structure formation dominated by cold dark matter 
(CDM), the current standard model of the cosmic struc- 
ture formation. The match of the power spectrum in the 
COBE scale with that in the galaxy clustering provides 
the most conspicuous evidence for the model. 

While power spectrum analyses are the most common 
way to demonstrate the consistency of the hierarchical 
clustering scenario, the amplitude estimated from galaxy 
clustering receives unknown biasing factors associated 
with galaxy formation mechanisms; therefore, this is not 
appropriate for a quantitative analysis as presented here. 
We consider matching of the two normalizations of the 
cosmic mass density fluctuation power, the normaliza- 
tions derived from COBE at a several hundred Mpc scale 
and the rich cluster abundance at z=0 which measures 
the power at « 8h~^ Mpc scale and most conveniently 
represented by the rms mass fluctuation parameter ag. 
We do not use the as parameters derived from velocity 
fields or other observations, which are more susceptible 
to various uncertainties. The advantage of using the clus- 
ter abundance information is that it refers to the mass 
function that is not affected by any biasing uncertain- 
ties, and the fiducial length scale 8 Mpc is close to 
that of clusters before collapse. The requirement of this 
matching leads us to derive quantitative constraints on 
the neutrino contribution to cosmic structure formation. 

We first consider the flat universe with low matter den- 
sity, but also discuss later the case for open universes. 
We assume a flat (Harrison-Zeldovich) initial power spec- 
trum, P{k) = Ak" with n — 1, which is the most natural 
prediction of inflation. The fluctuation spectrum receives 
a modification as Ak'^T(k) for a large k as the fluctua- 
tions evolve [16]. The shape of the transfer function T(k) 
depends on the assumed cosmological model, and the 
neutrino content. We use the computer code CMBFast 
[17] to calculate the transfer function for many choices of 
parameters. The spectrum is normalized with a fitting 
formula around £ = 10 deduced by Bunn & White [18] 
from the four-year COBE-DMR data [19]. They have es- 
timated one standard deviation error to be 7% in square 
root of the harmonic coefficient of CBR anisotropics Ce- 
We take the baryon fraction Hb = 0.015/i^ correspond- 
ing to 7710 = 4 [20]. The result is not very sensitive to 
this choice. 

We calculate the specific mass fluctuations within a 
sphere of a radius of 8 Mpc by integrating the spec- 
trum with a top hat window: 
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r<8h-^ 



dk4'Kk'^\6k\ 



sin(fcr) — (fcr) cos(A:r) 



• (1) 



neutrino mass density Q,^ is presented in Fig. 1 for a 
flat universe (a) Q = 0.3 and A = 0.7, and (b) ft = 
0.4 and A = 0.6. A set of curves (increasing towards 
the right) gives contours of constant (Tg. Another set of 
curves shows the neutrino mass density for given neutrino 
mass. 

The neutrino mass that concerns us is in the range < 
leV for most cases. The core radius of neutrino clustering 
allowed from the phase space argument [21] is 

= 3.2Mpc(m,VlcV)"2(i;/1000km s^^)^/^ (2) 

which is large compared with the core radius of rich clus- 
ters Rc ~ (0.12 ± 0.02)/i-iMpc-^ [22] for velocity dis- 
persion V « 10'^ km s""'^. Together with small neutrino 
mass density, we can ignore the neutrino component in 
integrating the cluster mass. We have estimated the con- 
tribution from neutrinos to the cluster mass within linear 
perturbation theory. The inclusion changes the result at 
most by a few percent, which can safely be ignored in the 
present argument. 
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FIG. 1. Contours of constant as derived from the COBE 
normalization in the Ho — ^Iv, plane for flat universes 

(O + A = 1). The region that satisfies the matching 
condition with the cluster abundance is indicated by a 
shade. Another set of curves indicates netrino mass density 
= (3mi,/93.84eV)/i-2. (a) O = 0.3, (b) = 0.4. 

This calculated as is compared with the value esti- 
mated from the rich cluster abundance. The estimate 
of as has been made by a number of authors [23-25,6]. 
The most ambiguous in such analyses is the estimate of 
the cluster mass, but the modern results are well con- 
verged among the authors, at least for 2; rs clusters. 
A summary is presented in Table 1. We take the values 
given by Eke et al. [24] which agree with other esti- 
mates within the error: as = 0.93 ± 0.07 for Cl = 0.3 and 
as = 0.80 ± 0.06 for Q. = 0.4. Adoption of Viana & Lid- 
die's [25] value makes the derived limit on neutrino mass 
slightly tighter. If we add the normalization error of the 
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CBR anisotropics in quadrature, the errors become 0.10 
and 0.09, respectively. The allowed range is shown by 
shadows in the figure. 

Wc sec from Fig. 1 that one can obtain the limit on 
neutrino mass if the Hubble constant is set. For Hq = 70 
km s-i Mpc-i we obtain 0.21 eV (r = O^/il < 5%) 
for VL = 0.3 and 0.91 eV (r < 15%) for Vl = 0.4. For 
Hq = 80 km s-i Mpc"^ our limits are 0.62 eV (r < 10%) 
for Q. = 0.3 and 1.8 eV (r < 22%) for O = 0.4.Our limit 
is summarized by a fitting formula: 

< [5.20/i(f2/0.3)2-03 - 3.20{n/0.3y-^^]h'^. (3) 

Allowing for conservative parameter space, fl < 0.4, 
Ho < 80 km s^i Mpc'^ and to > 11.5 Gyr, the up- 
per limit is 0.87 eV, which corresponds to r < 13% of the 
total mass density. 



(a) n = Q.3 open (b) = 0.4 open 
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FIG. 2. Same as Fig. 1, but for open universes. 

A similar figure is given in Fig. 2 for zero-A universes, 
(a) for n = 0.3 and (b) for fi = 0.4. The normaliza- 
tion from the cluster abundance is erg = 0.76 ± 0.09 and 
0.87±0.09 including the CBR normalization error. There 
is no consistent parameter range for — 0.3 for Hq < 100 
km s^^ Mpc^^ with or without neutrinos. A consistent 
parameter range appears for = 0.4, but only with a rel- 
atively high Ho- No- neutrino models are consistent for 
70 < Ho < 80 km s'^ Mpc'^. Requiring Ho < 80 leads 
to nil, < 0.5 eV, which is significantly stronger than the 
one for the flat universe. 

Tlic modification of power spectra with inclusion of 
massive neutrinos has been discussed by Hu et al. [26]. 
The change is about by a factor of two for rrii, = 1 
eV when Ho = 70 km s~^ Mpc~^. The scatter among 
the various data and our ignorance of the biasing factor 
make it difficult to exclude 1 eV neutrinos using the cur- 
rent power spectrum data. A strong constraint, such as 
nil, < 0.4 eV, would be obtained only when the power 
spectrum is derived from a homogeneous galaxy sam- 
ple with statistics as high as that would be expected in 



the Sloan Digital Sky Survey [26]. This would provide 
us with an alternative mean to set a limit on the neu- 
trino mass, although we must still assume the biasing 
factor being scale independent. The reason wc obtained 
a strong limit in this paper is ascribed to the advantage of 
using the cluster mass function, which is directly related 
to the mass fluctuation, as well as using the information 
spanning a very large baseline in the length scale. 

Let us discuss possible uncertainties or loop-holes of 
our argument. We have ignored the contribution from 
gravitational wave perturbations to CBR spectrum. Its 
inclusion only makes the limit on neutrino component 
more stringent. A possible loop-hole in our argument is 
the possibility that the index of the power spectrum is 
significantly larger than one. The COBE data alone do 
not exclude an index in the range 0.9 < n < 1.5, but 
the range is reduced to 1 < n < 1.2 if supplemented 
by other CBR data on small scales [12]. With an in- 
dex, n < 1, which can be easily realized with inflation 
models, the constraints become tighter. If n > 1, the 
excess large-scale power generated by neutrino perturba- 
tions arc cancelled by the intrinsically small large-scale 
power and more massive neutrinos become viable. For 
n = 1.2, the limit for Cl = 0.3 and h = 0.7 is loosened 
from 0.2 to 0.7 eV, and for h = 0.8, 0.6 to 1.4 eV. For 
our {Ho, n,to) range discussed above the hmit is 1.8 eV, 
still quite strong. We remark that we need some tricky 
tuning to give n > 1 in inflation models [27]. 

The limit we derived in this paper is quite strong. It 
is 5-20 times stronger than would be obtained from a 
straightforward mass density consideration < 93.8^lh^ 
eV. Ellis and Lola [2] have recently developed an argu- 
ment for neutrinos with a degenerate mass as large as 
5 eV with interesting physics. Such neutrinos, however, 
bring a large mismatch into the fluctuation power be- 
tween the very large scale and the cluster scale, causing 
a disaster to currently accepted cosmic structure forma- 
tion models. 

Let us finally compare our limits with those obtained 
from experiments or other cosmological considerations. 
A direct limit on electron neutrino mass from tritium 
beta decay is < 4.4eV (95% CL)[28] allowing for some 
systematic effects that make the measured negative. 
Additional limits are available if the neutrinos are of the 
Majorana type. The limit of lifetime for double beta de- 
cay of ^^Ge has now increased to 5.7 x 10^^ year (90% 
CL) [29] , which leads to 0.2 — 1.5 eV depending on the nu- 
clear matrix element used. We also refer to a limit from 
cosmological baryon excess: the condition for baryon 
asymmetry left-over leads to the Majorana neutrino to 
be < 1 — 2 eV [30]. The limit obtained in this paper 
does not depend on neutrino types. If one would con- 
sider only one species of neutrinos being massive, the 
mass limit simply becomes weaker by about a factor of 
three. 
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open 



1. 0.52 ± 0.04 0.80 ± 0.06 0.93 ± 0.07 0.76 ± 0.06 0.87 ± 0.07 

2. 0.56 0.86 0.99 0.77 0.84 

3. 0.53 0.79 0.80 
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